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We present a theoretical and numerical analysis of a quantum system that is capable of functioning 
as a heat engine. This system could be realized experimentally using cold bosonic atoms confined to 
a double well potential that is created by splitting a harmonic trap with a focused laser. The system 
shows thermalization, and can model a reversible heat engine cycle. This is the first demonstration 
of the operation of a heat engine with a finite quantum heat bath. 
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Recent developments in modeling the quantum dynam- 
ics of isolated systems [1] have provided a quantum un- 
derstanding of thermalization. However, they rely on ex- 
perimentally challenging systems. This has meant that 
development of a quantum understanding of heat engines 
has not been possible. In this Letter, we consider a model 
of cold bosonic atoms confined to a double well potential. 
Our analysis shows that this system exhibits thermaliza- 
tion when one well is initially more energetic than the 
other, and furthermore shows that the system can per- 
form a heat engine cycle. The system could be realized 
experimentally with present technology. 

The state of a quantum system evolves via the ap- 
pHcation of a unitary operator. Consequently, all pro- 
cesses are reversible and until recently it was not clear 
how a quantum system could reach thermal equilibrium. 
Srednicki suggested a solution called the eigenstate ther- 
malization hypothesis (ETH) Q , in which all eigenstates 
have the properties of a thermal state. Recently, this hy- 
pothesis was tested against other hypotheses by Rigol et 
al. [sl using hard core bosons confined to a lattice and 
found to be the only one that agreed. 

Several experimentally simpler schemes have been 
studied. Ponomarev et al. demonstrated thermal equili- 
bration when two isolated systems with different temper- 
atures were coupled through interactions between atoms 
in the different systems [4]. Chianca et al. allowed tun- 
neling of atoms between the systems in a four site lattice 
configuration [5] . However no simple measure of temper- 
ature could be obtained meaning it could not be devel- 
oped into a heat engine. A quantum description of a heat 
engine has been produced in Refs. [g] where the system 
is coupled to a classical heat bath. 

In this Letter, a quantum analysis of a heat engine 
is performed where both the system and heat bath are 
modeled quantum mechanically and atoms can tunnel 
between wells. Each well is described by two energy lev- 
els allowing temperature to be defined. One well acts as 
the system and the other serves as the heat bath. Heat 
is added and removed by tunneling of atoms between the 
wells. The volume of the wells are changed by modifying 
the harmonic confinement and the barrier position. All 




FIG. 1. Schematic of a double well created by splitting a 
harmonic potential with a focused laser. The diagram shows 
the possible tunneling and how the energy levels change due 
to the interactions. 



processes are done at a rate that provides good reversibil- 
ity of the engine. The system is modeled using both the 
truncated Wigner approach (TWA) to calculate the sys- 
tem for thousands of atoms and the accuracy is checked 
using calculations of the full quantum dynamics (FQD) 
for tens of atoms. Both results agree well showing an 
increase in number reduces fluctuations. 

N bosons confined to a double well trapping potential, 
Vdw(^), are described by the Hamiltonian 
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At low temperatures only the lowest laying single particle 
states are populated. Here we take into account the first 
two states as depicted in Fig.[TJ Therefore the field opera- 
tors can be described in terms of the four localized single 

particle functions, ip{x) = Yld^o + ^k(^)^k) ' 

where are the bosonic annihilation operators of an 
atom in well r and energy level / and described by the 
single particle functions This leads to the two-band 
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where we have ignored interactions between atoms in 
different wells. The ground and first excited state en- 
ergies are E[ = j dx^l^{x) (-I^V^ + VdU^)^ </>r(^)- 
The tunnel coupling between the wells is = 
J dx(t)^l{x) (-|^V^ H- Vdw{x)^ The interaction 

between atoms in the same well and on the same energy 
level is [/^ = g J and on different energy levels 

is U^^ = gj This term also leads to 

atoms changing energy levels. 

We consider a harmonic potential with oscillator fre- 
quency cjo, which is split by a focused laser beam located 
at xq from the center of the harmonic potential and de- 
scribed by a Gaussian potential Vq exp(— (x — xo)^/2cr^). 
The barrier height Vq = lOficJo, with width a = O.llho^ 
where Iho = y/h/mujo is the harmonic oscillator fre- 
quency. For a symmetric well, localized functions rep- 
resenting the energy levels in the different wells were cal- 
culated from the single particle eigenstates of the sys- 
tem. This gives J^/huJo = 0.153, J^/huJo = 0.226, 
E^/hwo = 1.37 and E^/hwo = 3.31. The interaction 
terms can be calculated from the integrals above and the 
interaction coupling, g. The interaction coupling can be 
varied by the Feshbach resonance [9] and for our purpose 
we use U^huJo = 2/7V, = SU^A and U^^ = /2. 
During the heat engine cycle only El. are modified, be- 
cause other terms are approximately constant. 

Thermalization of the system can be understood by 
first writing the initial wavefunction in terms of the 
eigenstates of the Hamiltonian, |^(0)) = 
where \(j)oi) are the eigenstates with energies Ecx and 
Ccx = (^al^(O)). The time evolution is given by, = 
e~*^^/^|^(0)). The expectation value of an observable O 
is thus given by 



id) 



a/3 



(3) 



where Oa^ = ((/)a|C^|(/>^). As stated above, ETH de- 
scribes the long time dynamics of the system [2]. The 
eigenstates of the system that make up the initial state 
have the properties of thermal state. After a consider- 
able time, the off-diagonal terms in Eq. (|3j) average to 
zero. This is called the long time average (LTA), which 
is a weighted average of the eigenstates and is indepen- 
dent of time, (O) = \Ca\'^Oaa- If the system is ther- 
malized, then the long time average should be equal to 
an ensemble average, which is given by a microcanonical 
average over a narrow energy interval around the mean 
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FIG. 2. The expectation values of the occupation numbers in 
an energy window near the mean energy of the system Eq for 
N — 40. The intersection with the mean energy of the sys- 
tem (dotted line) gives approximate values of the thermalized 
values. This supports the ETH. Inset: The dynamics of the 
field from t = to t = lOcj^^ in the symmetric double well 
for two slightly different initial conditions. The two trajecto- 
ries starting at A and ending at B(B') are not periodic, they 
gradually diverge and cover finite space. This demonstrates 
chaos in our system pj] ]. 



energy of the system, {d)Eo = ^I^Y.a,Eo±/^E/2 



Here M is the number of states in the energy interval 
A£^, £^0 = ^c^\Ca\^Eo,\s the mean energy of the system. 
We investigate the number of atoms in each of the energy 
levels. Therefore, the following equality is expected to be 
valid for thermalization 
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According to the ETH knowing a single eigenstate is suf- 
ficient to compute averages, O^a ^ {^)ec, [2, 3]. For 
O = nl this is shown in Fig. O The ETH also assumes 
initial conditions narrow in energy. For the initial state 
we use (see below) {{H'^) - {H)'^)^^'^ / {H) - 0.1. 

There are no general arguments supporting the ETH. 
There are proofs that quantum systems, whose classical 
analogues are chaotic, satisfy the ETH in the semiclassi- 
cal limit [10]. It is not computationally tractable to calcu- 
late the classical counterpart of our system for hundreds 
of atoms, therefore we look at the mean- field version of 
our system first. In this approach we replace the oper- 
ators with complex numbers b^ b^ in the Heisenberg 

equations of motion ^ = ^ [H^ b^] and get 
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where / = {1,2}, r = {L,R} and L = R, R 
dynamics of the system is chaotic (see e.g. 
example the imaginary and real parts of b\ shown in the 
inset of Fig. [2] clearly exhibit chaotic behavior. It follows 
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FIG. 3. (Color online) The heat engine cycle. TOP: the 
energy of the localized states in the double well system as 
the system evolves. BOTTOM: the population expectation 
value of each of the energy levels for TWA (dark lines) and 
FQD (faint lines). The initial state is |iV/2, 3iV/8, 0, 
and is first allowed to thermalize. The tunneling and in- 
teraction strengths are kept constant throughout the simu- 
lations at J^/hwQ = 0.153, J^/hwQ = 0.226, /fmjQ = 2/N, 
= 3U^/4 and [7°^ = 



that the trajectories are sensitive to initial conditions. 
We apply the truncated Wigner approximation, which is 
a powerful tool for the investigation of quantum dynam- 
ics in interacting many body systems ^3| . This amounts 
to solving Eq. (|5|) with different initial conditions, which 
accurately samples quantum noise in the system [l^ , and 
averaging afterwards. Since the trajectories of the sys- 
tem are chaotic, this might lead to equilibration within 
the TWA. For a coherent state, we sample the initial con- 
ditions as 6 = 6o + ^(i^i + 11^2) J where |6oP = n + l/2 
and Ui are Gaussian random variables. The coherent 
state corresponds to a condensed state of bosons. An- 
other possibility is a Fock state, which is sampled as 
b = {p -\- qv)e^'^^^ ^ where is a Gaussian random vari- 
able and ^ is a uniform random variable in the interval 

[0,1], P = 
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initial configuration of = 10^ atoms used is 

= 3A/'/8, = 0, = A'/S, so the right well is more 
energetic than the left well. The number of atoms in the 
excited state is small so the temperature is low enough to 
satisfy the two band approximation (see later). The sys- 
tem undergoes thermalization up to a time t = lOOco'^^ 
(see Fig. [3]). The expectation values of the occupation 
number of particles approach a constant value, which is 
independent of the initial state (being coherent state or 
Fock state) and agrees well with ETH. We find the sys- 
tem has thermalized with values nJ/A^ ^ 0.089 for both 



the left and right wells, and + = A^/2. 

Another method simulates the FQD by creating 
the Hamiltonian in the Fock basis of localized func- 
tions and propagates using the unitary matrix U{t) = 
ex.-p{—iHt/h). To make this approach computationally 
tractable the system is limited to 40 atoms. However, 
this is enough to verify the results of the TWA. Figure [3] 
shows that the FQD gives values that oscillate around the 
TWA results. These oscillations reduce for larger num- 
bers of atoms and for the N = 10^ atoms we expect good 
agreement with the TWA. Eq. (jlj is verified and com- 
pared with nj, which we average between times tuo = 80 
to tujQ = 100 to account for the oscillations, and gives 
good agreement, nl/N ^ Ifi^/N ^ {uDeo/N ^ 0.089. 

The heat engine scheme is similar to the Otto Cycle 
[lij . and was chosen to be experimentally simple. Heat 
is transferred from one well to another via tunneling of 
atoms. The work is done on the trapping potential by 
the expansion of atoms in a well when its size is slowly 
increased. Here the right well is the engine while the 
left well serves as the heat bath, and this corresponds 
to the clockwise solid line in Fig. HI The heat engine 
process starts at t = IOOuJq^. (a) {b): the laser barrier 
starts at xq = 0.2 /l^o and the harmonic confinement is 
adiabatically reduced to uj = 0.8cJo, so the populations 
of the localized functions remain roughly constant. This 
relates to the power stroke in the Otto cycle and work 
is done by the engine. (6) (c): the barrier moves to 
the left, Xq = — 0.2///io, leading to a rapid change in the 
populations. This is the heat rejection stage of the Otto 
cycle and heat is extracted from the right well, (c) {d): 
the harmonic confinement is adiabatically increased and 
the populations of the localized functions again remain 
roughly constant. This is the compression stage in the 
Otto cycle and work is done on the engine, (d) (a): the 
final stage sees the barrier move to the right and relates 
to the combustion stage of the Otto cycle. Here heat is 
added to the engine. While the volume does change in 
the steps (6) (c) and {d) (a), it is small relative 
to the adiabatic processes. The slight asymmetry of the 
evolution is due to the finite evolution time of the system, 
which was necessary to limit numerical inaccuracy. 

If a system is thermalized, it reaches a thermal state, 
which is described by the thermal density matrix pT = 
exp{-H/kBT)/TT{exp{-H/kBT)} [16]. The temper- 
ature T is chosen to best fit the data, and is used to 
calculate — Tr(pT^r)- At t = lOOcj^^ the system is 
thermalized with T\ ^ 2.QhjOQ/kB' We notice here that 
the two wells are not always in thermal equilibrium dur- 
ing the engine cycle, since we require the atoms flow from 
one well into another. However, the system reaches ther- 
mal equilibrium at t = 2006^^^ with T2 ~ 2.2h(jJo/kB and 
completes the cycle with T ^ Ti. The thermal energy, 
/c^T, is of the same order as the gap between the ground 
and first excited single particle states, — ^ 2Jiwo, 
which validates the two-level approximation. 
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FIG. 4. Pressure- volume plot for the right and left wells. Each 
well produces work proportional to the area enclosed by the 
curves. 



The pressure and volume (or length for a quasi-one 
dimensional system) in each well can be used to demon- 
strate the efficiency of the engine. If we assume the wells 
can be approximated by harmonic oscillators, then we 
can use the harmonic oscillator length as the length scale 
of the system, so Ir = ^JhjmuJr = h/ ^m{El — E^). 
For simplicity, we use the pressure for an ideal gas, 
P = 2Er/3V [17], where Er = n^^K + ^1^1 is the 
kinetic energy of a well, and taking V ^ 2/^, we get 
P = Er/3lr. The pressure and volume of the left and 
right wells are plotted in Fig. |4] and we see the hystere- 
sis of the heat engine, which shows work has been done. 
The whole isolated system cannot produce work, since 
the two processes in Fig. |4] cancel each other. This is 
consistent with the second law of thermodynamics. 

The entropy of a well is another quantity along with 
temperature, which can be measured to verify that the 
system has come into thermal equlibrium. It is given 
hy S = —k bTt {pr \og pr}^ where pr is the reduced den- 
sity matrix traced over the other well. As expected, it 
starts from zero at t = and approaches S ^ O.lksN 
at t = IOOuJq^. The entropy and other thermodynamic 
quantities can be measured from scanning the density 
profile of a quantum gas [18]. A focused electron beam 
with extremely high resolution of 150 nm can be used for 
this purpose [19]. 

The efficiency can be calculated as = W/ Qin, where 
work done is W ~ ^i=o i / '^r^^r heat added to 
the engine is Qin ^ ^i=o i / ^r^^r^ where = E^ ^ 
^ is the energy of a single particle shifted 
by the mean-field contribution from other particles in the 
same level [6]. We obtain 77 ^ 0.14. This is close to the 
Carnot efficiency r] ^ I - T2/T1 ^ 0.15. 

To experimentally realize the heat engine cycle, a mag- 
netic trap with radial and axial trapping frequencies of 
(jj ± / 27r = 4kHz and cJo/27r = 40Hz can be used to confine 
atoms in a ID trapping potential. In this case the ratio 
of the corresponding oscillator lengths is lho/l± ^ 10, 
such that the width of the focused laser beam is a ^ l±. 



For ^Li atoms this gives a ~ 0.6pm. A narrow laser 
beam was recently reported with a ^ 0.7pm and posi- 
tioning the beam to a lateral precision of 0.05pm 1201 . 
Taking the ID interaction strength g = 2h^as/ma\ [l7| 
yields the estimate for the two-body scattering length 
as ~ lOOnm/A/". The scattering length for ^Li atoms was 
achieved as small as ~ 10~^nm [21], therefore the num- 
ber of particles should not exceed N ~ 10^. For such a 
realistic setup the full cycle of the engine will take ~ Is. 

In conclusion, we have demonstrated thermalization of 
cold bosonic atoms confined to a double well potential, 
which is described by the two-band Hubbard model. We 
have shown that this can be utilized for realization of an 
isolated finite quantum engine, where one well produces 
work and another well serves as the heat bath. If realized 
experimentally, this will demonstrate the first isolated 
quantum engine with finite heat bath. 
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